Dynamical suppression of unwanted transition paths in multistate quantum systems 
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We introduce a method to suppress unwanted transition channels, even without knowing their 
couplings, and achieve perfect population transfer in multistate quantum systems by the applica- 
tion of composite pulse sequences. Unwanted transition paths may be present due to imperfect light 
polarization, stray electromagnetic fields, misalignment of quantization axis, spatial inhomogeneity 
of trapping fields, off-resonant couplings, etc. Compensation of simultaneous deviations in polar- 
ization, pulse area, and detuning is demonstrated. The accuracy, the flexibility and the robustness 
of this technique make it suitable for high-fidelity applications in quantum optics and quantum 
information processing. 

PACS numbers: 32.80.Qk, 32.80.Xx, 82.56.Jn, 42.50.Dv 



Introduction. Experiments in various branches of 
quantum physics require well-defined quantum energy 
states and well-defined interactions. For example, one 
of the basic conditions for the future quantum computer 
is a well defined qubit — a two-state quantum system. 
Real quantum systems, however, usually possess a huge 
number of quantum states and special care is needed to 
isolate just two of them. In real and artificial atoms 
this is usually done with polarized laser light, carefully 
aligned with the quantization axis. However, unwanted 
transition channels may still be present, which reduce the 
fidelity of the operations. For example, when an ultra- 
cold atomic ensemble held in a dipole trap is addressed 
by right circularly polarized (a + ) light, it is difficult to 
have all atoms "seeing" the same er + polarization be- 
cause not all of them are exactly in the focus of the laser 
fields; consequently, many atoms will "sec" an admix- 
ture of <7 + and o~ light (i.e., elliptically polarized light), 
which limits the fidelity of quantum state control. Un- 
wanted transition paths may be present also due to im- 
perfect polarization or alignment, stray electric and/or 
magnetic fields, off-resonant couplings, etc. 

In this paper, we propose a simple and efficient tech- 
nique for automatic compensation of such errors, which 
uses composite pulse sequences to dynamically suppress 
unwanted transition channels even without knowing the 
magnitude of these errors. We illustrate the technique in 
three-state and four-state quantum systems forming link- 
ages reminiscent of the letters "V" and "Y" , as shown in 
Fig. Q] the technique is, however, applicable also to more 
complex linkage patterns. We demonstrate compensa- 
tion of independent as well as simultaneous variations in 
polarization, pulse area and detuning. 

Composite pulse sequences have been used for several 
decades in nuclear magnetic resonance [H, 0] , and since 
recently, in quantum information processing and 
quantum optics as a versatile control tool for quan- 

tum systems. In fact, some of the basic ideas have been 
developed earlier in polarization optics in research on 
achromatic polarization retarders [ill EH ■ While com- 



FIG. 1: (color online) V system (left) and Y system (right). 



posite pulse sequences have been used mainly for two- 
state quantum systems, there are also studies of three- 
state and multistate systems 1(3, 13 1 . 

A composite pulse is a sequence of pulses with well 
defined relative phases, which are used as control param- 
eters in order to compensate imperfections in the excita- 
tion profile produced by a single pulse, or to shape the 
profile in a desired manner. The imperfections may be 
caused by an imprecise pulse area (due to fluctuating 
field intensity and pulse duration, spatial field inhomo- 
geneity, etc.), undesirable frequency offset (due to un- 
compensated electric and magnetic fields, dynamic Stark 
shifts, Dopplcr shifts, etc.), or unwanted frequency chirp. 

Here we use the toolbox of composite pulses to design 
recipes for compensation of imperfect light polarization, 
which may open unwanted transition paths and may turn 
a qubit or a simple three-state ladder into a complex tree 
of multiple states, with an unavoidable loss of efficiency. 

V system. The dynamics of a coherently driven V- 
shaped quantum system, shown in Fig. [1] (left), obeys the 
Schrodinger equation, 



\hd t c{t) = H(t)c(t), 



(1) 



where the vector (i) = [ci(t), C2(t), c^{t)\ T contains the 
probability amplitudes of the three states. The Hamilto- 
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nian in the rotating-wave approximation (RWA) reads 

H v (t) = (ft/2)A(n u - n 22 - n 33 ) 

+ (ft/2) [fiia We 1012 n 12 + n 13 (t)e^U la + h.c.] , 

(2) 

where A = loq — lu is the detuning between the laser 
carrier frequency uj and the Bohr transition frequency 
wo, and H^ = \j)(k\. The magnitudes of the Rabi fre- 
quencies are fljk(t) = \djk • E(i)|/ft, where E(f) is the 
envelope of the laser electric field and djk is the transi- 
tion dipole moment of the respective transition j <H> k: 
the phases of the Rabi frequencies are <f>i2 and <pi 3 . We 
assume that the two Rabi frequencies have the same time 
dependence f(t) and we introduce the root-mean-square 
(rms) peak Rabi frequency VL and the mixing angle via 
n 12 (t) = nf(t)cos9 and fl ls (t) = Qf(t)sm9. 

An important example of such a linkage is the transi- 
tion between the magnetic sublevel m = (state |1)) of 
a ground level with an angular momentum j = or 1, 
and the magnetic sublevels m = 1 (state |2)) and m = — 1 
(state 1 3)) of an excited level with an angular momentum 
j = 1 driven by an clliptically polarized laser pulse. The 
latter can be represented as a superposition of two circu- 
larly polarized a + and a~ pulses [14[: then, f2i2 = f2+ 
and f2i 3 = f2_, the angle of rotation of the polariza- 
tion ellipse is <f> = (4>i2 — </>i 3 )/2 and the ellipticity is 
e = cos 20. The values 6 = 0, tt/4, tt/2 (e = 1, 0, -1) cor- 
respond, respectively, to <r + , linear and o~ polarizations. 

Our objective is to transfer all population from state 
1) to state [2} and completely suppress the excitation 
channel |1) -o- |3). In the above example of magnetic 
sublevels, this can be achieved by a tr + polarized 7r pulse. 
However, if the polarization is not perfectly a + , then 
the atom will "see" some a~ polarized light and the un- 
wanted excitation channel |1) <H> |3) will open up. We 
show below that composite pulse sequences, instead of a 
single 7r pulse, can compensate such an admixture of un- 
wanted polarization, even without knowing its amount, 
and achieve perfect transfer efficiency. 

The V system described by the Hamiltonian (HI) ca n be 
transformed by the Morris- Shore transformation [15| into 
a decoupled state \d) = -e^ 13 sin0|2) + e" 1 ^ 12 cos0|3) 
and a two-state system composed of state |1) and a cou- 
pled state |c) = e 1 * 12 cos 6\2) + e 1 ^ 13 sin6»|3) driven by the 
following Hamiltonian 

H a (t) - (ft/2){A(n u - n cc ) + [0/(t)n lc + h.c.]}. (3) 

The corresponding propagator can be expressed in terms 
of the complex Cayley-Klein parameters a and b (with 
|a| 2 + |fe| 2 = l) as (3 



For resonant pulses (A = 0), with RMS area A — 
f* f fi/(i)t, the Cayley-Klein parameters are independent 
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FIG. 2: (color online) Transition probabilities Pi->2 (solid 
curves) and Pi->3 (dashed curves) for a V system vs the 
mixing angle 9 for a single pulse with RMS area it, and for 
composite sequences of three pulses (each with RMS area 7r) 
with phases </>i2 = (0,2/3, 0)n; <f>ri = (0, 1, l/3)7r, and five 
pulses with phases cj> 12 = (0, 1.411, 0.249, -0.432, -0.935)tt 
and 0i3 = (0, 0.454, -0.632, 0.14, -0.514)tt. 

of the pulse shape: a — cos A/2 and b = — i sin A/2. For 
A/0, a and b depend on the pulse shape. 

The propagator in the original basis reads [l(J 0, E3 

U(0) = 

a be^C be^S 

_ b * e -^i2 C a *C 2 + (S 2 (a* - e-*+SC , 
_ b * e -i<t>i3 S ( a *-C)e 2l *SC (C 2 +a*S 2 

(5) 

where S = sin6, C = cos 6 and C = exp[i f*f A(i)t/2]. 
Complete population transfer |1) — > |2) with a single 
pulse implies [ t/21 1 = lj which therefore requires a = 0, 
\b\ = 1 and 6 = 0. However, if 9 ^ 0, then the coupling 
between states |1) and |3) is nonzero and some popula- 
tion is unavoidably lost from state |2): either transferred 
to state 1 3) or left in state |1), because it is impossible to 
have \U 2 i\ = 1 for ^ (recall that |6| ^ 1). 

Deviation of from can be compensated to an ar- 
bitrary order by composite pulses. The propagator of a 
composite sequence of n pulses reads 

UW =U(0 n )---U(0 2 )U(0 1 ), (6) 

where 4>k = (4>i2 1 ^13^ ) are phase shifts of the fc-th pulse 
in the sequence with U(4>k) given by Eq. (f5]). The phases 
4>ij are free parameters. We choose to determine them 
by setting Pi->2 = \U^ '| 2 = 1 for = and nullifying 
the coefficients in the Taylor expansion of P\^,2 vs to 
the highest possible order. Since the global phase is ir- 
relevant, we take 4>i — without loss of generality. This 
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FIG. 3: (color online) Transition probability Pi_>2 for a V 
system. Upper frames: Pi^2 vs the mixing angle 6 and the 
RMS pulse area A for a single resonant pulse (upper left) and 
a composite sequence of three resonant pulses (upper right) 
with phases i2 = (0,2/3,0)vr; = (0, 1, 1/3)tv. Lower 
frames: Pn2 vs the mixing angle 6 and the single-photon 
detuning A for a single rectangular pulse of duration T and 
RMS area A — n (lower left) and a composite sequence of three 
rectangular pulses, each with duration T and RMS area tt, and 
phases <f>i2 = (0, 1/3, 0)7r; (f>i3 = (0, 2/3, 0)7r (lower right). 



compensation is demonstrated in Fig. [2] for composite se- 
quences of resonant pulses. For longer sequences (larger 
n), the transition profile Pi^2{0) widens and the un- 
wanted transition path |1) — > 1 3) is suppressed for a wider 
range of 9. Remarkably, for sufficiently long composite 
sequences, the transition |1) — > |3) can be suppressed 
even if its coupling is larger than that of the transition 
1) -> 1 2), i.e. in the range 9 > tt/4. 

Next, we have designed composite sequences which 
compensate simultaneous deviations in the RMS pulse 
area A from tt and the mixing angle 8 from 0. In order 
to find the composite phases we use the Taylor expan- 
sion of Pi_>.2 with respect to both A and 9, and annul 
the coefficients of as many successive terms as possible, 
while requiring also Pi_>.2 = 1 for A = tt and 9 = 0. This 
double compensation vs A and 9 is shown in Fig. [3] (top 
frames) for a sequence of three resonant pulses. 

In a similar manner, we have designed composite pulses 
which compensate simultaneous deviations in the mixing 
angle 9 from and the single-photon detuning A from 
resonance. This double compensation is demonstrated 
in Fig. [3] (bottom frames) for constant detuning: even a 
three-pulse composite sequence greatly expands the high- 
fidelity area around the desired point (9 = 0, A = 0). 

We note that even a triple compensation — vs 9, A and 
A — can be performed in a similar manner as above; it 



is, however, more difficult to illustrate it. 

Y system. The method for suppression of unwanted 
transition paths is readily extended to more complex sys- 
tems. Here we describe a natural extension of the V sys- 
tem to a Y-shaped system in which a state |0) linked to 
state |1) is added, as shown in Fig. Q] (right). The RWA 
Hamiltonian reads 



H y (i) = U v (t) - (h/2)AU 00 

+ (h/2) [n i(*)e'*' 1 IIoi+h.c.] 



(7) 



with Hy(f) being the Hamiltonian @ of the V system. 
The coupling of the additional transition |0) <r> |1) is pa- 
rameterized by a Rabi frequency with magnitude fioi(f)j 
which should share the same time dependence /(f) as 
the other two Rabi frequencies, and phase 0oi (which 
provides an additional control parameter). In addition 
to the mixing angle 9 in the V system, we introduce a 
second mixing angle £: f2oi(f) = ^ sm £/(f), ^12 (f) = 
f2 cos £ cos #/(£), and f2 13(f) = f2cos£sin#/(f), where 
now Qf(t) = ^fl 01 {t) 2 + n 12 (t) 2 + n 13 (t) 2 . Hereafter 
we take £ = tt/4, i.e., O i(f) 2 = ^i 2 (f) 2 + fMf) 2 - The 
couplings in the Y system in each interaction step can be 
caused by the simultaneous application of two lasers, one 
on the lower transition |0) — > |1) and another (elliptically 
polarized) on the upper V system |3) <— |1) — > |2). Al- 
ternatively, if in Fig. Q] we bring state |0) up to the set of 
states |2) and |3), the Y system can be viewed as a tripod 
system; if the latter is formed of the magnetic sublevels 
in the transition j = f> j = 1 then the three couplings 
can be produced by a pair of pulses, with elliptical and 
linear polarizations, derived from the same laser. 

The objective now is to transfer the population from 
state |0) to state |2) along the path |0) — > |1) |2), 
while suppressing the transition path |1) — > |3). Mathe- 
matically, this requires \U2o\ 2 = 1, Uqo = Uio = U^o = 0. 

As in the V system, because the couplings share the 
same time dependence /(f) and the Y system is on two- 
photon resonance, it can be transformed by a MS trans- 
formation into a set of two decoupled states and a two- 
state system. This allows us to obtain an exact ana- 
lytic expression for the propagator in the original basis 
in terms of the Cayley-Klein parameters of the MS two- 
state system, similar to the one of E q . © 
The propagator U for the fc-th pulse pair depends now 

on three phases: 0qi , 012 an< ^ 013 • 

Several conditions must be satisfied in order to achieve 
the desired transfer |0) — > |2). When the mixing angle is 
9 = (no unwanted coupling |1) — > |3)), this is achieved 
by a pair of simultaneous resonant pulses with RMS area 
A = J* f f2/(f)t = 27T As in the V system, nonzero 
9 can be compensated, even without knowing its value, 
by a composite sequence of pulses. 

The propagator of a sequence of n pulse pairs is given 
by Eq. ^ where now 4>k = {4>m \ ^12 \ ) ■ Compos- 
ite sequences, which compensate deviations of 9 and the 
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FIG. 4: (color online) Transition probabilities Po^2 
(solid curves) and Po->3 (dashed curves) in a Y sys- 
tem vs the mixing angle 8 for a single pulse pair with 
RMS area 27i, and for composite sequences of two and 
six pulse pairs (each with RMS area 7r), with phases: 
0oi = 012 = (0, 0) and 0i3 = (0, 1)tt for two 
pairs, and O1 = (0, 0, -0.181, -0.181, -0.033, -0.033)tt, 
0i2 = (0,0,-0.517,-0.517,-0.398,-0.398)77, and X3 = 
(0, 0.562, 0.026, -1.554, 0.393, 0.238)ti for six pairs. 



other interaction parameters, are constructed in the same 
manner as for the V system above: we expand the tran- 



sition probability Pq->2 
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in a Taylor series vs 



the relevant parameter (s) and determine the composite 
phases from the condition to annul as many terms (in 
ascending order) in this expansion as possible. 

Compensation vs the mixing angle 9 is demonstrated in 
Fig. H] for composite sequences of 2 and 6 pulse pairs. Si- 
multaneous compensation of deviations in both the mix- 
ing angle from and the RMS pulse area A from 7r is 
demonstrated in Fig. [5] (upper frames) for a composite 
sequence of 6 pulse pairs. Simultaneous compensation of 
deviations in the mixing angle 9 from and the single- 
photon detuning A from is demonstrated in Fig. [5] 
(lower frames) for a composite sequence of 6 pulse pairs. 

Conclusion. The proposed technique is a simple and 
efficient method for robust population transfer and sup- 
pression of unwanted transition channels in multistatc 
quantum systems. We have demonstrated this technique 
for three-state V and four-state Y systems, but it can 
readily be adapted to more complex systems. Unwanted 
transition channels may be merely unavoidable (e.g. due 
to off-resonant couplings), or can be activated, for in- 
stance, by deviations in light polarization or stray fields. 
By suitably choosing the phases of the constituent pulses, 
the unwanted transition channels are suppressed with 
very high fidelity even when the corresponding couplings 
are unknown, because the technique suppresses the un- 
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FIG. 5: (color online) Transition probability Po->2 for a Y 
system. Upper frames: Po^2 vs the mixing angle 6 and the 
RMS pulse area A for a single resonant pulse pair (upper left) 
and a composite sequence of six resonant pulse pairs (upper 
right) with phases O i = (0, 0, -0.986, -0.986, 0.348, 0.348)ti, 
0i2 = (0,0,0.667,0.667,-0.317,-0.317)71, and W = 
(0,-0.661,0.337,-0.042,0.955,0.285)tt. Lower frames: P ^ 2 
vs the mixing angle 9 and the constant single-photon detun- 
ing A for a single pair of rectangular pulses with RMS area 
A — 7i (lower left) and a composite sequence of 6 pairs 
of rectangular pulses, each with RMS area 7i, and phases 
0oi = (0,0,0,0,4/3,4/3)ti, 12 = (0,0,2/3,2/3,2/3,2/3)ti, 
and 0i3 = (0, 0.937, 0.854, 0.171, 0.798, 1.448)ti. 



wanted transitions for a broad range of coupling values. 
Compensation for deviations in the experimentally con- 
trollable parameters — polarizations, laser frequencies 
and laser intensities — can be done simultaneously with 
respect to several of them. The accuracy and the effi- 
ciency of the proposed technique and its experimental 
feasibility make it a potentially important tool in appli- 
cations requiring high fidelity of operation, such as quan- 
tum information processing and quantum optics. 
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